Charge Transport in the Dense Two-Dimensional Coulomb Gas 
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The dynamics of a globally neutral system of diffusing Coulomb charges in two dimensions, driven 
by an applied electric field, is studied in a wide temperature range around the Berezinskii-Kosterlitz- 
Thouless transition. I argue that the commonly accepted "free particle drift" mechanism of charge 
transport in this system is limited to relatively low particle densities. For higher densities, I propose 
a modified picture involving collective "partner transfer" between bound pairs. The new picture 
provides a natural explanation for recent experimental and numerical findings which deviate from 
standard theory. It also clarifies the origin of dynamical scaling in this context. 
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The Coulomb gas in two dimensions (2D CG) is the 
generic model for topological excitations ("vortices") in 
2D systems with a U(l) order parameter symmetry [jj]. 
Some important examples are the XY model of 2D planar 
magnets [||, superfiuid or superconducting films 
Josephson junction arrays B and, with some modifica- 
tions, 2D melting ElB. Vortices form a CG due to the 
characteristic, logarithmic form of their interaction po- 
tential. Furthermore, via a suitable choice of boundary 
conditions (an imposed current in the case of superfiuid 
systems) a topological "electric field" E can be applied 
on the vortices. For conceptual clarity, we shall use the 
CG language throughout this article ||. 

The globally neutral 2D CG in E = undergoes the 
famous Berezinskii-Kostcrlitz-Thoulcss (BKT) transition 
0. For temperatures T < T BKT , all particles are effec- 
tively bound in neutral pairs and the system is a dielectric 
insulator. At T BKT the pairs start to thermally dissoci- 
ate, leading to metallic screening and conduction. In 2D 
superfiuid systems, this unbinding of vortices causes the 
transition between superfiuid and normal phases. 

Besides its equilibrium phase transition, also the dy- 
namical behavior of the 2D CG is of considerable inter- 
est, since, e.g., in any superfiuid system, moving vortices 
are the main cause of dissipation. A few years after the 
BKT papers, a number of publications 1^-^] therefore 
extended the BKT pairing idea to nonequilibrium sit- 
uations, studying in particular the current response to 
a static, homogeneous applied field E. Their common 
procedure was to depart from the BKT decomposition 
of the system into two components, "bound pairs" and 
"free" particles, and to add kinetic equations for the par- 
ticle exchange between both; other interactions than the 
intra-pair one are cast into (both dielectric and metallic) 
screening of the latter. Ambegaokar et al. (AHNS) ||] 
extended the approach to an oscillating field E and con- 
cluded that at high driving frequencies the response is 
dominated by the internal dynamics of individual bound 
pairs (as well as, above T BKT , thermally generated free 
particles) whereas at low enough or zero frequency it is 
essentially due to free particle drift only. From BKT the- 
ory it is known |jj that below T BKT the charge correlation 



function decays algebraically as r P( T > with distance r 
where (3(T) := l/(e(T)T), T being the CG temperature 
and e(T) its static dielectric constant. The assumption 
of pure bound-pair reponse leads to algebraic decay in 
time as well with a dynamical exponent z — 2 at all 
temperatures T < T BKT ficfl . Field- induced pair disso- 
ciation allows for a nonvanishing static response below 
T BKT : AHNS found a nonlinear current- voltage (IV) re- 
lation j cx E a (J is the CG current density here H) with 



a(T)=p(T)/2 + l 



(1) 



Right at the BKT transition, /3(T BKT ) = 4 such that 
a(T~ KT ) = 3, which together with z = 2 satisfies the 
dynamical scaling ]Tl| ] relation a = z + 1. 

Away from T BKT however, doubt has recently been cast 
onto relatio n (|l| ) by dynamical simulations of the classical 
XY model ]12| as well as of the 2D CG in the continuous 
plane |l3| and on a lattice Q] , by a re-analysis of earlier 
experiments on layered superconductors Jf4[ , and also 
by new experiments on 2D proximity coupled Josephson 
junctions arrays [jl5|. A common conclusion of all of these 
investigations is that a rather seems to obey the relation 



a(T) = /3(T) - 1 , 



(2) 



proposed by Minnhagen et al. (MWJO) p2|. In other 
work [jl6|| better agreement with form ([!]) is found, so the 
issue is controversial. Right at T — T~ KT , both expres- 
sions give a value of 3, but (0) has a steeper increase with 
decreasing T. Assuming that the critical scaling relation 
a = z + 1 holds in an extended region below T BKT , MWJO 
link (0) to a temperature dependent dynamical exponent 



z(T) = P{T) 



(3) 



providing further scaling arguments in favour of (|^). 

In this letter, I propose a concrete physical picture for 
the correlated particle motion in a regime of intermediate 
CG particle density n, which in a natural way leads to 
both relations J9j independently. My picture is an 
extension of the "conventional" kinetic theory of Refs. 
|i|-[| and, accordingly, reproduces their results (|l[) and 
z = 2 for the "true" , asymptotic exponents in the limit 



1 



of small n. Furthermore, the underlying collective pro- 
cesses require a sufficiently strong applied field E (in the 
case of d.c. nonlinear response), or a sufficiently high fre- 
quency lj (in the case of a.c. linear response), respectively. 
Consequently, at fixed, finite n we also recover (Q) for 
small E, and z = 2 for small w. The detailed crossover 
conditions are derived below; their direct observation in 
experiments or simulations is however lacking so far. 

Following Refs. we consider overdamped particle 

motion, described by a system of Langevin equations 



dt 



ViVb (| aJ< - XjDqj 



E 



(4) 



where the CG particle positions and qi £ {±1} 

their charges (J^- qi = 0). fiQt Vo are the "bare" particle 
mobility and Coulomb potential; Vq behaves as Vo(r) ~ 
ln(r/ro) for r » rj, and Vq(t) ~ for r <C ro to prevent 
pair collapse jjj. The heat bath is represented by delta- 
correlated gaussian random forces £j. 

Let me briefly recall the conventional "free particle 
drift" picture J^-§] of IV response to a static applied 
field E. The dynamically screened pair problem is de- 
scribed by a two-particle Langevin equation of the form 
(||) with no replaced by an effective charge mobility fx e s, 
and Vo replaced by the statically screened pair potential 
V scr which we approximate as 



V (r)/e , r < 
V {l SCT )/e , r>£ 



scr 
scr 



(5) 



^ scr denotes the metallic screening length due to free 
charges and e the macroscopic, static dielectric constant 
of bound pairs up to this scale £ scr . Both /x e g and e are 
actually weakly scale dependent quantities, but this does 
not change the conclusions in any essential way and we 
ignore it here. h c h only sets the overall timescale and 
thus does not enter the exponents a or z. 

Let us assume T < T BKT (i.e., (3 > 4) for the moment, 
such that for E = all particles are bound in pairs and 
there is no metallic screening. Under the action of E ^ 0, 
pairs dissociate and free particles recombine with finite 
rates. Standard transition state theory yields for the re- 
spective rates per area Fdi ss cx yE° and r roc oc n 2 , where 
y is the fugacity of pairs at a distance r which controls 
the total particle density and rif is the density of free 
particles. In the stationary state, both rates are equal. 
Assuming then that the current is caused by free particle 
drift in the field E, we arrive at 



jdr oc n f E cx y 1/2 E p/ 



2+1 



(G) 



i.e., Odr = (3/2 + 1 [compare Eq. (EI)]. 

In order to see under which conditions the above ar- 
gument holds, let us look at length scales in the problem 
which are relevant to transport. One is the separation 
£e = 1/eE at which a pair becomes unstable in the ap- 
plied field E; pairs smaller than £e are (temporarily) 



bound. Another is the typical distance £<jr which a free 
particle drifts in the direction of E until it recombines, 



^dr := jdr/2r rec = jdr/2r d iss oc y 



-1/2^-/3/2+1 



(7) 



Since a neutral pair of size < £e is regarded as bound, 
a free particle density nf > £ E 2 or, equivalently, a "drift 
length" Ifc < £e clearly would be physically meaningless. 
Relations (g) and @) imply n f £ 2 E cx £ E /£dr oc y 1 / 2 E> 3 / 2 ~ 2 , 
so for sufficiently small y or E there is no contradiction. 
With increasing y and E, however, ml E and £e/ '£dr must 
exceed I at some point, implying a breakdown of the "free 
particle drift" picture and a crossover to a new "partner 
transfer" (PT) regime, where the constituents of a dis- 
sociated pair essentially have no opportunity anymore to 
drift freely, but immediately recombine with others. The 
PT regime may be visualized as a dense liquid of close- 
to-critical pairs of size ~ £e (by definition these can- 
not "overlap" each other) , "transferring partners" among 
each other by thermally activated hopping |l7| . The 
whole is of course immersed in a screening medium of 
smaller pairs, as discussed above. 

In order to estimate the current in the PT regime, let 
us focus on the typical distance £ which a given parti- 
cle separating from its partner is dragged along by E 
until it binds to the next. This length shrinks from 
^ ~ ^dr 3> (e in the "drift" regime to its lower bound 
£ rT ~ x-p T £e (with a dimensionless constant x rT > I) in 
the PT regime. Assuming then that pair dissociation in 
the PT regime is still approximately described by transi- 
tion state theory (implying Tdiss oc yE 13 ), we find 



Jv 



= 2F 



diss^PT y 



(8) 



i.e 



£Jpt = P — 1 precisely as obtained in Refs. |l2j |fj] 
[see Eq. (Q)]! Note that an extension of the Kostcrlitz 
renormalization procedure by including £ as a scaling 
variable |8p^| does not have a comparable effect, which 
indicates that the PT picture captures additional aspects 
of the strongly correlated dynamics at higher density. 

MWJO have attempted to reconcile the new rela- 
tion a PT = f3 — 1 with the "free particle drift" picture, 
postulating a modified dependence r roc cx n\ +2 ^^ ^ of 
the recombination rate on the free particle density rif. 
If however, as we have argued here, the concept of free 
particles itself becomes questionable in the new regime, 
this type of explanation clearly cannot be adequate. 

In the high temperature phase (J3 < 4) thermal pair 
dissociation leads to metallic screening, which provides a 
further characteristic length Q 



(9) 



As long as £ scr > £e , field induced pair dissociation dom- 
inates over thermal one and our arguments remain essen- 
tially valid, only the order of the "drift" and PT regimes 
is reversed: the PT regime now is at small E (but still 
intermediate y). At still smaller E such that £ scx < £e 
one eventually obtains an ohmic regime. 
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The above results for j(y,E) can be expressed by a 
scaling form of the nonlinear CG conductivity er CG 0, 

a CG (y, E) := j(y, E)/E = £ E 2 d(£ E /£ d ^ i E /£ SCI ) , (10) 

where the scaling function a depends on (3 (but not y) 
and has the asymptotic forms cx(w, 0) oc u for u <C 1 
("drift") and oc u 2 for u 3> 1 (PT), respectively, and 
ct(u, v) oc v 2 for any u and t> ^> 1 (ohmic regime). 

Estimating quantitatively the prefactors in Eqs. (0) 
and (|), we find that the new PT regime extends to rel- 
atively small y and E, so its predominance in experi- 
ments and computer simulations is not surprising. Using 
transition state theory to consistently estimate both pair 
dissociation P] and recombination H rates we obtain 



(11) 



As to the screening length £ SCI , a simple Debye-Hiickel 
estimate jl| turns out to be quantitatively dissatisfactory 

— 1/2 

since it implies £ SCI < n f ' at p « 4. Minimizing the 
grand canonical potential of a mixed system of bound 
pairs and free particles with respect to ^ scr yields instead 




-2/(4-/3) 



(12) 



Equations (]1 1|) , ( Jl 2| ) imply an extension of the PT regime 
to arbitrarily large £e (i.e., small E) in the vicinity of the 
BKT transition, since £ scr diverges as (3 — ► 4 _ and, for 
&ny > e~ 4 and any finite £e, £dr/£E shrinks to below 
1 as (3 —> 4+. For typical 2D vortex systems |l|Q|§) 
y and (3 obey to a good approximation the "Villain" re- 
lation In t/ w -/37r/2 ||. Figure @ shows the resulting 
location of the different regimes in the ft-E~ x plane for 
some reasonable values of the (so far unknown) constant 
particular, in the vicinity of the BKT transition 
the asymptotic exponent = (3/2 + 1 appears only at 
extremely small driving fields E which probably are not 
reached in these investigations. 

We now turn to dynamical correlations at E = 0. Fol- 
lowing Refs. we assumed here that the motion of 
a given bound pair or free particle is essentially uncorre- 
cted with the rest of the system, feeling it only through 
its screening properties. In particular, at T < T BKT , 
charge density correlations are governed by the internal 
dynamics of pairs JTo| . The effective dynamical exponent 
z for fluctuations with wavenumbers around k can then 
be read off from the typical relaxation time r(r) oc r z of 
pairs of size roc fc _1 ; it thus will depend on whether the 
relaxation process is dominated by "drift" or PT. The 
two mechanisms are illustrated in Fig. |^: (a) a member 
of a pair can simply drift towards the other under the ac- 
tion of its field, or (b) each of them can recombine with a 
member of another, smaller pair sitting in between them. 



In case (a) each particle has to move a distance oc r pulled 
by the field Ep oc 1/r of its partner, which implies 



Tdr(r) oc r/Ep oc r 2 



(13) 



2. This relaxation mechanism has been an- 



alyzed in detail by Ambegaokar and Teitel |10| and is 
certainly the dominant one in the limit of large r. In 
case (b) there is an area oc r 2 within the big pair in 
which its dipole field has a strength Ep, oc 1/r and thus 
induces the dissociation of smaller pairs with a rate per 
area Tdiss yE^, so the typical waiting time until such 
an event occurs behaves as Ml 



? (r) oc l/r 2 r diss oc y^r 13 " 2 



(14) 



i.e., z PT = (3 — 2 in agreement with (g). The four-body 
correlations implicit in this argument again underline the 
collective character of the PT mechanism. Comparing 
both relaxation times, Td r (r)/r PT (r) oc yr 4- ' 3 , shows that 
below the BKT transition (fi > 4) the "drift" mecha- 
nism is dominant at small y or k (oc r -1 ), whereas with 
increasing y and k, PT takes over at some point. 

( |l3| ) , (|lj) imply dynamical scaling forms of the charge 
density correlation function which in turn, via the fluctu- 
ation-dissipation theorem, suggest a dynamical extension 
of ( [To| ) to a CG (y, E,uj), including a further scaling vari- 
able £dr/£uj (£lu '■= \J HckT/uj). For sufficiently large y 
and u>, the PT mechanism dominates, leading to a sim- 



plified form cr CG (y, E, u>) 



'(£E ZpT /£i) (the form of 



the scaling variable reflects the fact that r PT (£e) oc £ dr . oc 
£~ ZrT in the PT regime). The validity of the "critical" 
scaling relation a P 

T — P T 

+ 1 in a whole T region below 
T BKT (which was an assumption in the MWJO scaling 
argument!) now follows very naturally by the usual rea- 
soning |y]]. It expresses the fact that in the PT regime 
both dc charge transport and the decay of charge fluctua- 
tions are governed by the same mechanism, characterized 
by a single relevant length scale £e oc 1/E. 

In summary, I have described a new dynamical regime 
of the 2D CG in which charge transport predominantly 
occurs through direct particle transfer among bound 
pairs, instead of free-particle drift. It offers a natural in- 
terpretation of recent experimental and numerical find- 
ings on 2D vortex dynamics which deviate from "stan- 
dard" AHNS theory @. In contrast to the more ex- 
treme view held by MWJO I do, however, recover 
AHNS behavior asymptotically for low driving fields and 
frequencies; the new picture thus resolves the apparent 
contradiction between the MWJO and AHNS scenarios. 

Further tests of this picture in experiments and nu- 
merical simulations are clearly needed; in particular, a 
direct observation of the predicted crossovers would be 
desirable. Another testable signature is the changed de- 
pendence of dc response (||) and relaxation times ([li]) on 
the pair fugacity y. I am not aware of any controlled way 
to tune y independently of (3 in real superconducting sys- 
tems. In experiments on superfluid 4 He films, the vortex 
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fugacity can be increased above the "Villain" value by 
adding 3 He impurities, which stick to the vortex cores 
and lower their energy pl| . However, as we have seen, 
in order to observe the PT/ "drift" crossover a reduction 
of the fugacity would be needed instead. 

In any event, in real systems a reliable determination 
of the effective CG temperature T and thus of usually is 
an extremely difficult task. In contrast, CG simulations 
are natural test candidates since here T, y are separately 
controllable model parameters (within certain limits this 
is also true for generalized XY type models p2f). Pre- 
liminary Monte Carlo data for a lattice CG p3] indeed 
support the proportionality (||) between j and y in the 
accessible parameter range. 

The crossover at low E [|| in the IV characteristics 
should in principle be observable both in experiments 
and simulations. Figure |l| suggests to look for it at tem- 
peratures sufficiently far below T BKT (/? > 5, say), since 
the crossover value of Ie diverges at T BKT . The same is 
true in principle for the low u> crossover in the dynamical 
correlations (i.e., in the ac linear response) although this 
one may be covered by other collective effects ]22| , |24| ]. 
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FIG. 1. Different IV regimes in the /3-_E _1 plane ac- 
cording to Eqs. ([nj), (^) and using the "Villain" relation 
In y — —f3n/2, for some reasonable values of x PT . The PT 
regime, characterized by the condition 4r/^PT < (.e < (.sa, 
expands with increasing x PT (shading is for x PT = 2). 




r r/2 
FIG. 2. Pair equilibration by (a) "drift" and (b) "partner 
transfer" (PT). In (b), the relevant region of small-pair disso- 
ciation, leading to PT relaxation of the big one, is framed by 
a dashed line. Electric field lines of the big pair are shown. 
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